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Abstract 

The direct scheme of Model Reference Adaptive Control 
System (MRACS) may have several disadvantages; (i) the 
current state of the plant cannot be easily grasped from 
values of adjustable adaptive controller parameters, (ii) it 
has often inferior conditions to an indirect scheme, with 
respect to the persistent excitation (PE) conditions of 
adjustable controller parameters for convergence to optimal 
values. Indirect MRACS based on Dynamic Certainty 
Equivalent (DyCE) principle hasnot been well studied 
though the importance of indirect scheme is recognized in 
an actual field using adaptive control. This paper presents 
new two design schemes of indirect MRACS based on DyCE 
principle. These are realized by recursive identifiers and new 
high order tuners. 
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Introduction 

The design method of Model Reference Adaptive 
Control System (MRACS) based on the Surrogate 
Model Control (SMC) scheme is originally proposed 
by Morse (Morse, 1992). The SMC obtained by a 
control strategy that the output of a certain identifier is 
always equivalent to the desired value as an output of 
reference model of control system, has several 
advantages as follows. It has the ability to apply 
MRAC to a certain plant; for example, a plant with 
nonlinearities which do not satisfy globally Lipchitz 
conditions, a LTI plant whose relative degree is n*or 
n* + 1, n* being a known positive integer and so on. 
Additionally, MRACS based on SMC scheme is 
superior in a transient response of controlled variable 



to the conventional MRACS based on certainty 
equivalence (CE) principle since the SMC input 
includes an auxiliary input which prevents the 
derivative of adjustable controller parameter from 
appearing in the tracking error. Hence, the swapping 
term in the tracking error in MRACS with CE control 
is perfectly removed by the auxiliary input. The 
adaptive control input by SMC scheme can be 
expressed as a result in a form that the auxiliary input 
is added to the CE adaptive control input. Therefore, it 
is also called an adaptive control based on the 
dynamic certainty equivalence (DyCE) principle. 
However, it needs high order differential values of 
adjustable parameters which must be updated by an 
adaptation law. Morse (Morse, 1992) and Ortega 
(Ortega, 1993) proposed a series of the high order 
tuner (HOT) being special adaptation laws, which can 
generate high order differential values of adjustable 
controller parameters though it does not use high 
order differential values of the control variable. 

Several modified schemes of MRACS due to DyCE 
principle have been proposed [for example, (Itamiya, 
1998), (Masuda, 2000), (Tanahashi, 2007), (Tanahashi, 
2011)] in consideration of the realization of HOT, some 
improvement of control performance and so on. 
However, the scheme based on SMC in which 
estimates of plant parameters are used directly for 
adaptive controller hasnot been well studied. In 
conventional MRACS based on CE control, such a 
scheme is well known as indirect one. Advantages in 
indirect method are (i) one can observe or predict a 
behaviour of plant through estimates, (ii) fast 
convergence property of estimates and small 
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parameter error may be expected since PE degree for 
an adaptation law is less than that in a direct MRACS. 
Especially, the advantage (i) can is never obtained in 
any direct MRACS but it is significant from practical 
viewpoints. 

Two new schemes of MRACS with DyCE control have 
been proposed here using estimates of plant 
parameters directly. These can be constructed by 
introducing certain recursive identifiers and HOTs 
induced as steepest descent methods which minimize 
certain integral costs. These HOTs are similar to the 
adaptive law proposed by Kreisselmeier 
(Kreisselmeier, 1977). In the schemes proposed here, 
any design calculation with nonlinear algebra isnot 
needed. The controllers are constructed directly by 
estimates of plant parameters. 

This paper is organized as follows; the next section 
states the system representation and the problem 
statement. The third Section presents the scheme 1 to 
design a MRACS with DyCE control directly using 
estimates of plant parameters and its brief stability 
analysis. The scheme 2 such that the steady state 
response of control variable is improved and its 
stability analysis is shown in the fourth section. Simple 
numerical simulation results in order to verify the 
effectiveness of the proposed schemes are indicated in 
the fifth section. The last section concludes the paper. 

Notations: The following notations are applied; s 
denotes Laplace operator. The symbol T expresses the 
transpose of a vector or a matrix. A definition symbol 
is expressed as: =and IHI means Euqulidean norm or 2 
norm of Matrix induced from Euclidean norm. The 
symbol L 2 an d L m mean the set of square integrable 
functions and uniformly bounded functions, 
respectively. Furthermore, 

(G(s)[v])(t):= f 0(t-T)t7(T)dT 



xW(£): = 



((s + Ay'M)(t)for; < 



[(D + X) ] x(t) for j > 

where G(s) is the Laplace transforms of g(t), jmeans 
integer, A is a positive constant and D is a time 

differentialoperator ^ and =^jx(t) norm is 
defined as 

II x t \\ 2S -.= ^e- s ^ \\x(t) II 2 dr) 

for x(-) £ L 2e and 8 > 0. L 2e means the set of square 



integrable function in finite interval. || W (s) H^, is the 
H^norm of a proper and stable transfer matrix W(s). 
In addition, || W(s) || 2 is the H 2 norm of a strictly 
proper and stable transfer matrix W(s). 

System Representation and Problem 
Statement 

The control object (plant) considered here is a LTI, 
SISO and minimum phase system. It has a coprime 
transfer function N ? (s)/D ? (s) where degrees of 
polynomials D p (s) and JV p (s) are n and n-n* 
respectively. When u{t) is the plant input in time t, the 
plant output y(t) can be represented as 

y(t) = T ((t) + 6(t) (1) 

where e(t) is a exponentially decaying term whose 
size depends on the initial value of the plant initial 
state, the plant parameter 8 £ R 2n_n * +1 and C(0 e 
E 2n-n*+i are defined by 



(2) 



9'- = [Su> Sy> SyV 

8 U : = [8 n », 8 n » +1 , 8 n ] T 

8 y : = [8 n +l> 9n+2> ■■■ > $n+n*-l] T ' 
8y\ = [8 n+n *, 8 n+n * + 1 , ... , #2n] T - 

[Cu T (0- £(0. CyW 

( u (t):= [u [ " n * ] (0. u[-( n * +1 M(t) ul-"] (t)] 1 



W-= \y^Kt),y^Kt) y l ~ {n - 1)J (t)]' 

r- = [yl- n *](t), yK n * +1 )l(t) y [_n] (t)] T 



(3) 



Remark 1: For any X > , D p (s) and iV p (s) can be 
generally represented as 

D p (s)=I l ? =0 a i (s+A)"-; a = 1(4) 

JVpO) = Hn-bi (s+A)"-'(5) 

Then, the next relationship holds; 



8 i = b i 
8 n+i = -a.{ 



(i = n* ~ n) 



(6) 



(i = 1 ~ n) 
The followings are assumed for the plant; 

Assumption 1: 

Al) The system degree n and the relative degree n* are 

known a priori. 
AT) 8 is unknown. 

A3) The sign of 8 n * is known a priori. Here, 8 n * > 

without loss of generality. 
A4) Available signals are only u(t) and y(t)- 

The problem to be considered is to design the MRACS 
with DyCE control which uses estimates of plant 
parameters. The desired output that the plant output 
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y(£) must follow is given by the reference model 
output y m (t); 



*.(«:- (g§[r])W<7) 



where r(t) is a piece- wise continuous and bounded 
reference input, D m (s) is a monic and Hurwitz 
polynomial with degree n m , iV m (s) is a polynomial 
whose degree is less than or equal to n m — n* . The 
transfer function is the desired reference model for the 
closed loop system. 

Scheme 1 

Design ofMRACS 

Now, consider the following linear recursive identifier 
corresponding to (1); 

yil (t):=0 T (tX(O ) 
y iy (t): = 9\t)xPj(t) (j = 2 ~ n* - 1) 1 (8) 

yin*(t):=£ T (t)</v(t) J 

where 0(t) £ R 2n " n * +1 is the estimate of 9, ^y(t) Q = 
2 ~ n*) is defined as follows; 

ty(0= = [£(0. &(0. &(0. (0] T ) 

Cyy(t):= [y [ - 1] (t),y h2] (t) yt-^Ct)] 7 [(9) 

^n*(t):=[^(t),Cj„*(0,(J(t)] T ) 

Cyn*(t):= [yi„*ii.yi„*- 2 yii V) 

Since yiy(t) is j times differentiable on t if (fc) (t) (fc = 
~ j) is available, here an adaptive control input can 
be proposed so as to satisfy the SMC law; 

5W(0 = ttn(0(H) 

The control law (11) leads to 

y(t) = y(t)-y m (t) = ^: i £iy(t)(i2) 

where £iy(t) is defined as 

£n(t): = y(t)-yn(t) 

eijCt): = 9ij-i(t) - 9ij(t) (j = 2~ n*)}^> 
The equation (11) means 

r'(t)-0 r (t)^r ] (t)-A(t) 



u(t) = 



B .(t) 



(14) 



where 0(t) 6 R 2n_n * and xj){t) £ R 2n_n * are sub-vectors 
which consist of elements exception of the first 
component in 0(t) and ip n *(t) respectively. r'(t) and 
/i (t) are defined as 



my = zjU „. Cy ■ (0«(t)) T i/#*~ ; ' ] (O 



(15) 



Remark 2: It can be shown that (14) is equivalent to the 
adaptive law used in a usual indirect MRACS if 

flit) EE 0. 

Remark 3: In a direct scheme MRACS with SMC, the 
control input can be realized because the regressor 
vector is n* times differentiable. However, the 
surrogate model control 

5>ii(0 = y m (t)(i6) 

cannot be realized since £^(t) (j = 2 ~ n*) cannot be 
obtained without using y(t) ~ y (n * _1) (t). 

The control law (14) needs the estimate (t) of plant 
parameter and its high order differential values. 
Therefore, we aim at the following relationship; 

£ll (t) = T (tX(t) + e(t) ] 

£ij(t) = ii:ie n m)4?-i(t)y 

where 0(t): =9 - 0(t) and j = 2- n*. From (17) and 
(12), it is clear that £jy(t) for j = 2 ~ n* converges to 
zero and then also y(t) converges to zero if n+ ;(-) £ 
Loo for i = 1 ~ n* — 1 and £n(t) vanishes. The 
following HOT can guarantee such a property; 

[HOT1] 

d(t) = r 1 (t){qi(t)-fi 1 (t)}(18) 
where 

r!(t): = r 10 P!(t)(i9) 

r 10 : = diag[y ln .,diag{y ln . +1 ,y ln . +2 , ...,y 12 „}](20) 
P(t): = diag{£ n .(t),/ 2n _ n .}/W 10 (t)(21) 



N 10 (ty.= ffi(t) + 2n-n*(22) 
<?i( t ):=(7^[yN^])(0(23) 
fi i(0:=(j^?[<N<S])(0(24) 



where A x > and yij > (j = n* ~ 2n) . yn(t): = 
y(t)/W 1 (t) and Cn(0 : = CCO/^iCO mean normalized 
signals where iV x (t) is defined as follows; 

Wi(t): = VPi + C T (tX(0; Pi > 0(25) 

Initial estimates are given as < 0„*(O) < oo and 
II 0(0) ||< oo. Note that from (18)-(24), high order 
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differential values of 0(0 can be obtained without 
using high order differential values of y(t). 

Remark4: It can be seen that 

*(o-«i(o*(o = K- Sf) T(26) 

hit): = Jo K{t - T){y N (r) - T (O C N (T)} 2 dr(27) 
where fti(t) means the inverse Laplace transform of 

Remark 5: The term, {q x (0 - R t (t) 0(0 } T , is the 
steepest descent direction of A(t) for 0(0 . It is 
resemble to one in Kreisselmeier's adaptive law 
(Kreisselmeier, 1977). The definition of h(t) > 
guarantees the existence of high order differential 
values of 0(0- 

Furthermore, n *(O in the definition (19) of T(t) 
contributes to the fact n * (t) > for any non-negative 
t (Sawada, 2013). Hence, conventional (switching) 
projection algorithm[(Goodwin,1987), (Ikhouane,1998), 
(Nagurney, 1996), (Kuhnen, 2004)] is not needed to 
satisfy n «(O > 0. 

The normalizing signals N 10 (t) and iV x (t) establish 
II i\(t) ||< 1 and |y N (t)| <ll || +c, II ( N (t) ||< 1. These 
results contribute to P^-), q^-), R^-) E L m . Thus, 
0(0 £ and the upper bound of II 0(0 || can be 
designed by T 10 if - 0(-) £ Loo- 
After all, the schemel proposed here consists of the 
control law (14) and the HOT (18). 

Stability ofMRACS Based on Scheme 1 

Firstly, the stability of the adaptive loop has been 
shown by schemel. The next lemma holds. 

Lemma 1 (Stability of adaptive loop): The HOT (18) of 
scheme 1 satisfies the following properties; 
Pl-1) 6 Lm _ 

Pl-2)A(-)eii- 
Pl-3) §(-)EL 2 . 
Pl-4)0 n «(O > Oforallt> 0. 

Pl-5) £ 11N (-)ei 2 

where £ 11N (Q: = £n(0/JVi(Q. 

proof: Let the positive definite function 1^(0(0) be 



Of. = 0t - 0;(O(28) 



which is defined in the set; 

C:= {0(0 £ R 2n - n * +1 |<9 n .(t) > 0}(29) 

Then, the time differentiation of Vi(0(O) evaluated 
along the trajectory of (18) becomes 

VdHW = -AWMoW < 0(30) 

Therefore, ^(0(0), 0(0 G Loo and ^(0(0) converges 
some constant in [0, Vi(0(O))]. Hence, Pl-1) is proven 
because q 0) (-), ^(0 G Loo for = ~ n*. From the 
integration of (30), Pl-2) is proven. Pl-3) holds since 
II 0(0 || 2 can be evaluated as 

II 0(t) ll 2 <ll I\ || 2 {trace ff 1 (0} 2 A(t)(31) 
and /,(■) £ L x . Pl-4) is shown by the fact 
lim§ n , ( tHo "co ^(0(^(0 " «i(0 6(0} = 0(32) 

where n c0 : = [-1, T ] T is the normal vector at n *(O = 
of the convex set C. Also, Pl-5) is proven by Pl-2) 
and the mathematical induction similar to the paper 
(Tanahashi, 2011). Q. E. D. 

From the above lemma 1, the following theorem about 
the main control loop (closed loop) holds. 

Theorem 1 (Stability): All variables in the adaptive 
control system which consists of (1), (14) and (18) are 
uniformly bounded in t, and 

lim^oo y(t) = 0(33) 

proof: In the following, the exponential decaying term 
is abbreviated because it does not affect the stability. 
Note that £(0 can be represented as 

at) = (2r(s)[y])(0 
7 rcv = r Dp(s) Dp(s) Dp(s) 

1 J - W) n *' (s+A)"* +1 ' "" (s+A)"' >(34) 



1 1 



(s+A) nJ 



s+A' (s+A) 2 

which satisfies (12) and z(s) can be decomposed into 

z(s) = b + z(s) 
b:= [1/0 



+ 1(35) 

n*> ®2n-n'] > 



where z(s) is a stable and strictly proper transfer 
vector. 

Then, the following relationship can be obtained. 

n* 

M(0 at) = (z(5)[y m ])(0 + b £ £ iy (0 

;'=2 
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n* 




+ U0) 


7 =1 


j (t)(36) 



where M(t) is the upper triangular matrix as 
/? T (t) 



M(t): = 



2n-n*^2n-n* 



(37) 



/?(t): = [0 n .(O- -0 T (t)] T /6n*(38) 

which is clearly non-singular from Pl-4) of the lemma 
1. Therefore, 

II C(0 ll 2 < c + cTf =2 1 £iy(t)l 2 +clf =1 II s ljt || 2 a (39) 

where c > is a generic constant. In the derivation of 
above inequality, the next relationships are used; 

II (z(s)[£ iy ])(t) || 2 <|| z(s) || 2 5 -|| e 1Jt ll 2 5 (41) 

where II z(s) || 25 means the H 2 norm of z(s — 5/2). 5 
satisfies < S < 8 ; 8 : = min{2/l, 2A Np }. A Np means 
mirij{ | Re. (5 Np j)|} where 5 Npi (i = 1 ~ n*) is a zero of 
iV p (s). 

Moreover, from (17) and §(■) £ L^, e i; -(t) (j = 2 ~ n*) 
is evaluated as 

K<t)\ 2 <C II £ iy _ lt || 2 5 | 

II £i/t ll| fi < c II £ iy - lt ll| J ( 

where c > is a generic constant, by using the 
concept of H 2S norm and H^g norm (Ioannou, 1996). 
Therefore, there exists in some generic constant c > 
such that 



\£ij(t)\ 2 <c II £ llt \\ 2 2S 
II £iyt l&< c II e llt \\ 2 2S 



Inequalities (39), (43), the definition of JV^t) and the 
property Pl-5) lead to 

II ftt) H 2 < C + C J C £l 2 1N (T) || C(T) || 2 dT(44) 

This is equivalent to 

II C(0 II 2 <c exp{c (T)dr}(45) 

from Bellman-Gronwall lemma (Ioannou, 1996). This 
means CO) £ L OT from Pl-5). Hence, y(-) £ L^, from (1). 
u(-) £ L OT is guaranteed from linear boundedness 
condition (Narendra, 1989) or regularity (Sastry, 1989) 
since the controlled object is a minimum phase system, 
y(-) £ Loo and the fact that u(t) grows at most 
exponentially from <?(•) £ Loo- 



The result of u(-), y(-) £ Loo means £ 1X (0 £ n L 2 
and <f(0 £ Loo- So that, £n(t) is uniformly continuous 

signal from 0(0,^(0, CO), CO) e ^ and 

£ll (t) = T (t)C(t)+£ T (t)C(t)(46) 

These facts lead to zero convergence of £n(t) by 
Barbalat's lemma. Therefore, y(t) converges to zero 
from (17) and (12). Q. E. D. 

Scheme 2 

The scheme 1 proposed in the previous section is able 
to achieve the property that MRACS is stable and the 
tracking error y(t) onverges to zero. However, y(t) 
and £n (t) can never go to zero simultaneously. 
Meanwhile, y(t) may have ill convergence property if 
6 n+ i(t) for i = 1 ~ n* — 1 remain as large values since 
the tracking error equation is described by (12) with 
(17). In order to reveal such a phenomenon, another 
scheme is proposed here. 

Design of MRACS 

Consider another recursive identifiers corresponding 
to (1) as 

y 21 (t):=0 T (t)C(O 

y 2 j(t) = y 2 y-i(t) 

-«y-i(t){y [ - (y - 1)] (0 - £ T (t) C [ -°'- 1)] (t)} \(47) 

ay-i(t): = XtX+i (0 «y-i-i(0 ; « (0 = 1 

Q = 2 ~ n*)> 

Then, y 2 y(t) in (47) is j times continuously 
differentiable without the derivatives of y(t) . For 
example, when n = n* = 2, y 22 (t) is evaluated from 
(47) as follows; 

y 22 (t) = S 2 (t)u[- 2 ](t) 

+e 3 (t)£ T (t)C hl] (t) + e 4 (t)y [ - 2] (0(48) 

It is clearly 2 times continuously differentiable without 
derivatives of y(t) if 0(t), 0(t) and 0(t) are available. 

As the adaptive control law in scheme, the next SMC 
law has been put forward; 

ym<t) = y m (t)(49) 

The control input is obtained by solving the following 
equation on u(t); 

in- ] (t) = y l f\t)(50) 

According to (49), the tracking error y(t): = y(t) - 
y m (t) becomes 
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m = m-o T (t)cp(t)(5i) 

= T (t)<KO + e(t)(52) 
where 

m-- = z£?«t (0y H] (0(53) 

*»(*): = Z?=^«( (t)C H] (0(54) 
eW-S&^i (t)e H] (t)(55) 

Since y(t) in (52) is affine on 9(t), and e(t) decays 
exponentially if #„+;(■) (i = 1 ~ n* — 1) £ L M , it means 
that y(t) goes to zero if 9 T (t) cp(t) + e(t) tends to zero 
as t increases. Therefore, 9(t) and its high order 
differential values can be updated by 

[HOT 2] 

S(O = r 2 (t){q 2 (t)-S 2 (t)0(t)}(56) 
where 

r 2 (t):=r 20 P 2 (t)(57) 

r 20 : = diag[y 2n ., diag{y 2n . +1 , y 2n « +2 , ... , y 22 „}](58) 
P 2 (t): = diag{0 n .(t),/ 2n _ n .}/W 2O (t)(59) 

JV 20 (t): = J^.(t) + 2n-n*(60) 
^(0:= (^F[y N <P N ]) (0(61) 

fl2 W := (^F [<PN ^ ] ) (t)(62) 

where A 2 > and y 2 j > (j = n* ~ 2ri) are adaptive 
gains. <Pn(t):= (p(t)/N 2 (t) means normalized signals 
where N 2 (0 is defined as follows; 

jv 2 (t): = Vp 2 + <p t (0<p(0; P2 > 0(63) 

Initial estimates are given as < n «(O) < oo, || 0(0) ||< 

CO. 

Note that from (56)-(63), high order differential values 
of 9(t) can be obtained without using high order 
differential values of y(0- 

Clearly, high order differential values of 9(t) can be 
obtained from (56) without derivatives of y(0- 

Remark 6: (56) is induced from the steepest descent 
method with smooth projection to minimize the 
following cost function; 

M0 = Jo h 2 {t - T){y N (r) - 9\t) «p N (r)} 2 dr(64) 

where h 2 (t) means the inverse Laplace transform of 
Af/(s + A 2 )"*. 



The total surrogate model controller using plant 
parameters estimates based on the scheme 2 consists 
of the control law (49) and the HOT (56). 

Stability ofMRACS Based on Scheme 2 

Stability of MRACS with the scheme 2 is shown as 
follows. (56) satisfies the next lemma. 

Lemma 2 (Stability of adaptive loop): The HOT (56) of 
scheme 2 satisfies the following properties; 

P2-1) £ L m . 

P2-2)/ 2 (-) £ Ll 

P2-3) £(•) £L 2 . 

P2-4) <9 n «(t) > Oforallt> 0. 

P2-5) y N (0 e L 2 where y N (0 = = Kt)/N 2 (t). 



proof: The proof can be shown in the same way as one 
in the lemma 1. Let the positive definite function 
V 2 0(t)) be 

2n 

+ 2y7^ (t); ^ : =^"^(0(65) 

i=n*+l 2i 

which is defined in the same set as (29). Then, the time 
differentiation of V 2 (9(t) evaluated along the 
trajectory of (56) becomes 

^(0(0) = -; 2 (o/w 2 o (o< o(66) 

Therefore, V 2 (9(-)), £(■) £ L m and V 2 (fl(0) converge 
some constant in [0, V 2 (8(0))]. Hence, P2-1) is proven 
because q 2 \-)> R 2 \-) £ L°o for j = ~ n* . From the 
integration of (66), P2-2) is proven. P2-3) holds since 
II 9(t) || 2 can be evaluated as 

II 0(0 H 2 <ll r 2 (0 II 2 {trace R 2 (t)fJ 2 (t)(67) 
and / 2 (-) £ l x . P2-4) is shown by the fact 
lim§ n , (tH0 "co r 2 (0{<? 2 (0 " ^(00(0} = 0(68) 

where n c0 : = [-1, T ] T is the normal vector at 9 n *(t) = 
of the convex set C. In addition, P2-5) is proven by 
P2-2) and the mathematical induction similar to the 
paper (Tanahashi, 2011). Q. E. D. 

From the above lemma 2, the following theorem about 
the main control loop (closed loop) holds. 



Theorem 2 (Stability): All variables in the adaptive 
control system which consists of (1), (49) and (56) are 
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uniformly bounded in t, and 
lim^„ y(t) = 0(69) 



proof: In the following, the exponential decaying term 
is abbreviated. Note that (p(t) can be rewritten as 

<Kt)=I?=oV (0(w t (*)[y])(0) 

, x i , x (70) 
W;(S): = tz(s) 

where z(s) is the same transfer vector as (34). 

By (52) and procedure similar to the previous section, 
<p{t) can be represented as follows; 

M(t)cp(t) = c m (t) + O(s)[0> N ■ JV 2 ])(t) 

n*-l 

+ £ ffi (t)(w i (s)[0 T <p N ■ JV 2 ])(t) (71) 

i=0 

n*-l 

c m (t): = ^ «j (t)(w t (s)[yU)(t) (72) 



i=0 

Where M(t) and z(s) are defined as (37) and (35) 
respectively. 

Then, || cp(t) || 2 is evaluated as 

II <p(t) || 2 < c + c jjyg (t) II <p(r) II 2 dr(73) 

Since M -1 (-) £ L^, «;(■) £ L m , and P2-5) are satisfied. 
c > is some generic constant. This means 

||<p(t) || 2 <c exp{c f y^(x)dr}(7A) 

from Bellman-Gronwall lemma (Ioannou, 1996). This 
means cp(-) £ L m from P2-5). Hence, y(-) £ from (1). 
u(-) £ Loo is guaranteed from linear boundedness 
condition (Narendra, 1989) or regularity (Sastry, 1989) 
since the controlled object is a minimum phase system, 
y(-) £ Loo and the fact that u(t) grows at most 
exponentially from §(■) e L m . 

The result of u(-), y(-) £ Loo means y(-) £ Loo n L 2 and 
<p(-) £ Loo- So that, y(t) is uniformly continuous signal 
from §(;),§(>), (p (•).?»(•) 6 Loo and 

y(t) = ^ T (t)«p(0 + e T (t)«p(t)(75) 

By Barbalat's lemma, these facts lead to zero 
convergence of y(t)- Q- E. D. 

Numerical Example 

In this section, simple examples will be presented to 
illustrate the usefulness of the theoretical results of the 
proposed two schemes. 



The followings are settings for the numerical examples. 
The plant was chosen as 



[u] (t) 



0.5s 

and the reference model was selected as 
1 



y m (t) = 



s 2 + 1.4s + 1 



[r] (0 



where r(t) is the rectangular wave with the period 
20[s], varying between and 1. Other settings are as 
follows; 

A 1 = 1,A 1 =A 2 = 2,p 1 = p 2 = 0.1, 

r 10 = r 20 = 207 3 ,e(0) = [2, o, of 

Then, the plant parameter 8 = [1, 2.5, — 1.5] T . 

Fig.l and Fig. 2 are simulation results. These illustrate 
the usefulness of the theoretical results of the 
proposed schemes. It can be seen that the tracking 
error in MRACS with the scheme 2 is small relative to 
one with the scheme 1. 



'a -2 




s 

5 



t[ S ] 




FIG. 1 SIMULATION RESULTS BASED ON SCHEME 1 
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Conclusions 

New two schemes for an indirect MRACS based on 
SMC are proposed. It is indicated that recursive 
identifiers and HOTs with smooth projection function 
for # n *(t) > in t > play central roles in both 
scheme. Extensions to the robust version to 
unmodeled dynamics and disturbances are the focus 
of our future studies. 




+», 0.4 ■ 
« 0.3 ■ 

0.2 - 



20 40 €B SO t [s] 

FIG. 2 SIMULATION RESULTS BASED ON SCHEME 2 
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